Spectral fluctuation properties of constrained 

unitary ensembles of Gaussian— distributed random 

d^ : matrices 
o 
o _ 

^ , Z Pluhaf ^ and H A Weidenmiiller^ 

Ch ■ 1 

. Institute of Particle and Nuclear Physics, Faculty of Mathematics and Physics, 

Charles University, 18000 Praha 8, Czech Republic 
OO ■ 2 Max-Planck-Institut fiir Kernphysik, Postfach 103980, 69029 Heidelberg, Germany 



> 
00 

cn 



E-mail: Hans . WeidenmuellerOmpi-hd . mpg . de 



Qh' Abstract. We investigate the spectral fluctuation properties of constrained 

(-H I ensembles of random matrices (defined by the condition that a number Nq of matrix 

elements vanish identically; that condition is imposed in unitarily invariant form) in 
the limit of large matrix dimension. We show that as long as Nq is smaller than a 
critical value (at which the quadratic level repulsion of the Gaussian unitary ensemble 
of random matrices may be destroyed) all spectral fluctuation measures have the same 
form as for the Gaussian unitary ensemble. 
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■ 1. Introduction 

^. 

5t , We investigate the spectral fluctuation properties of the constrained unitary ensembles 
of Gaussian-distributed random matrices (CGUE) introduced in Ref. p. Constrained 
ensembles of random matrices deserve interest because they represent entire classes of 
non-canonical random-matrix ensembles that were proposed to mimic typical properties 
of interacting many-fermion systems more closely than do the canonical ensembles 
(the Gaussian Orthogonal, Unitary, and Symplectic Ensembles) |2]. Examples for 
constrained ensembles are the Embedded Gaussian Orthogonal Ensemble [H H] and 
the Two-Body Random Ensemble [5[ El [7] . The constraints essentially require certain 
matrix elements or linear combinations of matrix elements to vanish. It is difficult 
to deal with the constrained ensembles analytically because they lack the invariance 
properties that give analytical access to the canonical ensembles. That difficulty is 
overcome by imposing the condition of unitary, orthogonal, or symplectic invariance 
on the constrained ensembles. It is expected that that condition leaves the spectral 
fluctuations unchanged (whereas the eigenfunctions acquire the same distribution as for 
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the canonical ensembles). Here we focus on the case of unitary invariance, i.e., on the 
CGUE. 

Some spectral properties of the CGUE were exhibited in Ref . [I] . In particular, the 
following sufficient condition for level repulsion was given. The quadratic level repulsion 
characteristic of the CUE (the ensemble of Gaussian-distributed unitary random 
matrices) also prevails for the CGUE provided that the number Nq of constraints does 
not exceed a critical value, 

Nq < AT-* (1) 

with Nq^^ defined in Eq. (fT6l) below. 

In the present paper we go beyond Ref. pj and address the spectral fiuctuation 
properties of the CGUE in the limit of large matrix dimension ^ 1. We do so for 
Nq < N^'\ For Nq = the CGUE coincides with the CUE. For Nq > N^'^ the 
form of the constraints does not seem to permit definitive analytical statements. It has 
remained an open question to what extent the spectral fiuctuation properties of the 
CGUE (beyond the statement of sheer level repulsion) are the same or differ from those 
of the CUE for < A^q < Nq'^. We prove the following 

Theorem: For matrix dimension ^ 1 and the number of constraints Nq < Ag''* 
(with A'q'* defined in Eq. ( fT6l) below), the spectral fiuctuation measures of the CGUE 
coincide with those of the GUE save for correction terms of order 1/A^. 

To make the paper self-contained, we collect in Section [2] some definitions and 
results given in Ref. [I]. In Section [3] we slightly modify the definition of the constrained 
ensembles so as to remove a singularity. In Section H] we discuss the form of the 
constraints in the limit A^ ^ 1. Our proof is given in Section [5l It is based on an 
approach developed in Ref. [8]. Section [6] contains a discussion. Some technical details 
are presented in an Appendix. 



2. Definitions 



We consider Hermitean matrices acting on a Hilbert space Ti. of dimension A^. For any 
two such matrices A and B, we introduce the canonical scalar product in terms of the 
trace 



{A\B) = Ti{AB) . (2) 

us 

which obey 



This allows us to define an orthonormal basis of A^^ Hermitean basis matrices Ba = -B^ 



{B^\Bp) = Ti{B^Bp) = (3) 

and 

Y,\B^){B^\ = 1n (4) 

a=l 
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where Iat is the unit matrix in N dimensions. Any Hermitcan matrix H acting on Ti. 
can be expanded in terms of the N'^ Hermitean basis matrices as 

H=J2 haBa. (5) 
a=l 

The Gaussian Unitary Ensemble (GUE) of random matrices is obtained by assuming 
that the expansion coefficients ha are uncorrelated Gaussian-distributed real random 
variables with mean value zero and a common variance. For the GUE, the probability 
density W(H) of the matrix elements of H has the form 

iy(//)d[//] = A^GUEg^p f _ JL(^H\H))d[H] . (6) 



Here d[H] is the product of the differentials of all independent matrix elements, jV^^UE 
is a normalization factor, and 2A is the radius of Wigner's semicircle. 

We introduce the constraints by considering two orthogonal subspaces labeled V 
and Q with dimensions Np and Nq — N'^ — Np, respectively. These are defined in terms 
of orthogonal projection operators 

Np 

V= Y.\Bp){B,\ , 

Q= E \B,){B,\. (7) 



q=Np+l 



We have 



Constraints can be formulated in the form {H\Q) — for all H. In the CGUE such 
constraints are used in unitarily invariant form. The CGUE is defined by the probability 
density W-p of the matrix elements of H given by 

Wr{H)d[H] =A^GUEg^p ^_^^H\H)y[H] 



X /d[C/](n ^i\l^2iUB,U^\H))) . (9) 

The integral d[U] extends over the unitary group in dimensions. The Haar measure 
of the unitary group is normalized to one, i.e.. 



d[U] = 1 . (10) 
We diagonalize the matrix H with the help of a unitary matrix V, 

H = VxV^ , (11) 

where x = diag(a;i, . . . ,xn) is the diagonal matrix of the eigenvalues. The integration 
measure becomes 

d[H] oc A\x)d[x]d[V] , (12) 
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where dx is the product of the differentials of the eigenvalues, where d[V^] is the 
Haar measure of the unitary group in N dimensions, and where A(a;) denotes the 
Vandermonde determinant 

A(x) = n (^A* - • (13) 

l<lM<i/<N 

Eq. f|T2l) shows that eigenvalues and eigenvectors of the CGUE are uncorrelated random 
variables. The joint probability distribution P-p(x) of the eigenvalues is given by 

Pp(x) = AToexp (^-^{x\x)^A\x)Fv{H) , (14) 

where 

FAH)^ lm{m,l^^{B,\UHU^))'^ 

is a function of the eigenvalues {a;^} only, and where A/q is another irrelevant 
normalization factor. Comparison of Eq. ([9]) and Eq. ([6]) shows that the eigenvalue 
distribution of the CGUE differs from that of the CUE by the factor F-p(if). GUE-type 
level repulsion is contained in the factor A'^{x) in Eq. (IT^ . and such level repulsion will 
prevail also in the CGUE unless F'p{H) is singular whenever two eigenvalues coincide. 
In Ref. [1] it was shown that F-pi^H) cannot be singular if the number Nq of constraints 
obeys the inequality 

Nq < AT-* = N{N - l)/2 - X: %(% - l)/2 • (16) 

i=i 

Here it is assumed that the matrix B = J2q SqBg with real coefficients Sg possesses 
asymptotically (all Sq large) J sets of degenerate eigenvalues with multiplicities Lj, 
j = l,...,J. 

3. Modified Form of the Constraints 

The function F'p{H) embodies the constraints. Therefore, it is the central object of 
study in this paper. The treatment of F-p{H) simplifies when all confining matrices Bq 
are traceless. We believe that that case is physically the more interesting one, for the 
following reason. We show in the Appendix that whenever the BqS are not traceless, 
there always exists an orthogonal transformation of the set {Bg} to a new set {Bq} 
such that F-p{H) is unchanged and that all Bq with q > Np + 1 are traceless. The one 
constraining matrix Bj^p^i that is not traceless, is the sum of a traceless part and of 
a multiple of the unit matrix. But constraining H with a multiple of the unit matrix 
means that we constrain the centroid of the spectrum of H. We cannot think of a 
physically interesting situation where such a constraint would be meaningful. This is 
why we focus attention on the case where all Bq are traceless, 

{Bq) = for all g = A^p + 1, A^P + 2, . . . , A^^ ^ (17) 
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and treat the more general case where the condition fll7l) is violated, in the Appendix. 
Here and in the sequel we use the symbol (A) to denote the trace of the matrix A. This 
is consistent with the definition ([2]). 

For the developments in Section [5] we note the following properties of F-p{H). 
The function F'p{H) is real (this follows from Eq. (ITSi) ) and positive definite (this is 
seen when we write the defining delta functions as limits of Gaussians). Using Fourier 
transformation, we can write F-p{H) in Eq. f|T5|) as an iVg-fold Fourier integral, 

9=1 ■ 



/ \ 2 X Nq/2 ^^Q i- . 

^^^^^ = (2^) n J d5,y d[U] exp{^{Bis)\UHU^)} (18) 



where 



Bis) = J2s,B, (19) 



F (m^( V""^' /d. detexp{zx^&,(g)} 
^^(^^ °^ l2^j J A(x)A(b(s)) ■ 



and where s stands for the set {si, . . . , siy^}. The integral over the unitary group can 
be worked out and with ds = Jig dsg yields (see Ref . P ) 

^ detexp{ix^6^(s)} 
A{x)A{b{s)) 

Here the b,^{s) are the eigenvalues of the matrix B{s), and A{b) is the Vandermonde 
determinant of the b,y{s), see Eq. (1131) . In Ref. [Ij it was shown that the integrals over 
s converge if the condition (fT6|) is met. Moreover, inspection of Eq. (l20l) shows that 
F-p{H) is not singular when two eigenvalues x^,Xu coincide. 

However, because of the form of the constraints (Eq. f|T5l) ). the function F-pi^H) is 
singular when all eigenvalues of H coincide. To see this we define 

H = H-'^{H), (21) 

use the assumption f|T71) and rewrite Eq. f|T8l) in the form 

Ft.{H) = Fr{H) = ^^Y""^' J ds J d[U] exp{z{B{s)\UHU^)} . (22) 

Eq. fl22]) shows that F-p{H) is singular when H = 0, i.e., when all eigenvalues of H 
coincide. The singularity mirrors a singularity in the definition (fTSjl of Fp{H). Indeed, 
when x^ = x^ = ?/ for all /i, = 1, . . . , A^, the Dirac deltas in the definition (IT^ take 
the form 6{y{Bq)). Because of Eq. f|T7|) each of these terms is singular. We avoid the 
singularity by modifying the definition of F'p{H). Instead of F'p{H), we consider the 
constraining function 

(^2)^iVQ/2 

iVA2 

The factor in front of F-p{H) guarantees that F-p[H) is not singular at if = 0. At 
the same time, that factor is a function of the sum of the eigenvalues only. Thus, 
that factor cannot modify the correlations of close-lying eigenvalues x^, and the spectral 
fiuctuation properties of the constrained ensembles defined by the constraining functions 
F-p^H) and F'p{H) are the same. Moreover, for real the function F'p{H) is real and 



Fv{H)=(\t^) F-p{H). (23) 
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positive definite. According to Eq. fl20p Fp{H) is not singular for finite values of the 
Xfj^. Inspection shows that when one of the eigenvalues, say, tends to infinity, F-p{H) 
cannot grow more strongly than some power of x^. That growth is much weaker than 
the Gaussian suppression of large eigenvalues in Eq. (Q. Hence the confinement of 
the spectrum to a finite interval characteristic of the GUE persists also for the CGUE 
with constraining function Fp{H) although the shape of the average spectrum may be 
modified. 

Collecting everything, we have 

^^^^^ = {^^V I "^[^1 / ^^^{'^UB{s)U^\H)] . (24) 
It is convenient to introduce the new variables tq = Xsq. Then 

Fv{H) = {^^Y'" j m J dt eM^{UB{t)U^\{H/X))} . (25) 

This shows that F-p{H) depends on H only via the dimensionless ratio H j A, as expected. 
Because of unitary invariance, F'p{H) can depend only on unitary invariants constructed 
from H/\. The only such invariants are the normalized traces of if" /A" with positive 
integer n. For ^ 1 this is shown explicitly in the next Section. The probability 
density for the Hamiltonian matrices of the CGUE is given by 

Wv{H)d[H] = ATexp ^- — {H\H)yr{H) d[H] . (26) 

The substitution of F-p{H) by F'p{H) also modifies the normalization factor of W-p but 
that is irrelevant for what follows. 



4. Asymptotic Form of F-p{H) for A^ > 1 

For A^ ^ 1 we now display explicitly the dependence of F'p{H) on the normalized 
unitary invariants {1/N){H^/X) with positive integer n. We mention in passing that 
the terms of leading order in a systematic expansion of Fp{H) in powers of Nq/N"^ can 
also be obtained from the Harish-Chandra Itzykson Zuber integral 0,110], or from the 
standard supersymmetry approach [HI [12]. We use the assumption (fT7|) and discuss the 
case where not all constraining matrices are traceless in the Appendix. Then {B{t)) = 0. 
We consider the expressions 

/ d[U] {z{{H/X)\UB{t)U^)r (27) 

with k positive integer. Such expressions are generated when the exponential in Eq. (l25l) 
is expanded in a Taylor series. To calculate the integral over the unitary group we use 
a method valid for A^ ^ 1 [T3l HI]. To leading order in 1/A^ the integral can be done 
using Wick contraction on the matrices U, the rules being 

U^Ml il/N)S^J,, and U^Mpa . (28) 

Terms of higher order are obtained in a similar fashion and lead to similar results 
but are not considered here. We first look at a few simple cases. For k = 1 the 
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expression (I27|) vanishes. For k = 2 and = 3 we obtain (i'^ /N)[{1/N){{H/ X)"^)] 
{B^{t)) and 2\{i^ /NY[{1/N){{H/ Xy^)]{B^{t)) , respectively. For A; = 4 Wick contraction 
generates two terms. One is proportional to the square of the k = 2 term just considered. 
The other is given by 3\{i^/N^)[{l/N){{H/\y)]{B^{t)). For k = 5, Wick contraction 
generates two types of terms: The product of the k = 2 term and the k = 3 term, and 
a new term given by 4! (iViV^) [{1/N){{H / X)^)] {B^{t)). 

For the general expression fl27j) we consider all partitions of k into sets of 
positive integers ki, k2, . . . kf greater than unity such that J2i=iki = k. To leading 
order in expression (1271) is given by the sum over all such partitions, the 

contribution of each partition being i'' Q^) {^l^') x ... x (^^''''"^■'^f-''^ lli=i{ki - 
l)\{l/N''^-^)[{l/N){{H/X)''^)]{B''^{t)). The terms of higher order in also involve 
products of traces of powers of H/X and of traces of B, the difference being that at 
least one trace of a power of H/X is multiplied by at least two traces of powers of B 
such that the sum of the exponents of B equals the exponent of H/ X. It is shown below 
that (-B") and nj(-B"*) with ^iTii = n are of the same order in N so the neglect of such 
terms is legitimate. 

We conclude that to leading order in the integral over the unitary group in 
Eq. fl25|) is given by 

d[U] e^v{^{UB{t)U^{H/X))} 

= exp{^(l/n)(^7iV«-i)[(l/iV)((i7/A)-)](i?"(t))} . (29) 

n>2 

For ^ 1 the constraining function F is then given by 

^^^^^ - [2^0^) 

X fdt exp{Y^{l/n){^yN-''m/N){{H/Xr)]{B^^m . (30) 

n>2 

It may seem that because of the factors A^"~^ the terms of higher order in n in 
Eqs. (!29|) and (l30l) can be neglected. We now show that for Nq ~ A^^ this is not the 
case. In Eq. ( l30l) we replace the Cartesian integration variables tq by polar coordinates 
{r, fi} in Nq dimensions where 

r^ = E^? (31) 

and where Q stands for the angular variables. We write 

B{t) = rB{n) (32) 
and since {Bq\Bqi) = 6qqi have 

{B'm = 1 . (33) 

Let b^{fl) denote the real eigenvalues of B{Q). Then I]^&^(f2) = 1 and |&;x(f2)| < 1 
for all = 1, . . . , A^. For integer n > 2 this implies that 

(5"(1])) < 1. (34) 
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This, incidentally, justifies the omission of terms of order above and shows that 
{B"'{Q)) and ((if/A)") are characteristically different: The first expression is (at most) 
of order unity while the second is of order A^. That is why we always carry the second 
expression in the form {1 / N) {{H / X)"") . 

Using the transformation to polar coordinates we observe that the term with n = 2 
in Eq. (|30|) gives a Gaussian integral in r. Expanding the terms with > 2 in the 
exponent in a Taylor series we are led to consider radial integrals of the form 

J dr r^«-i+2'=exp{-cr2} . (35) 

Here c is a constant and 2k must be even as otherwise the integrals vanish. Compared to 
the leading term {k = 0) these integrals yield a factor (A^q+2/c— 2) {NQ+2k—4) x . . . x Nq. 
(We assume for simplicity that Nq is even). If the expansion of the exponential converges 
sufficiently rapidly so that for Nq ~ N"^ we need consider only terms with k <^ Nq then 
every power of r in the exponential in Eq. fl30l) effectively carries a factor ^JNq, and 
the series in n proceeds effectively in powers of Nq/N"^. For Nq ~ A^^ that factor is of 
order unity. 

While it is, thus, not permitted to terminate for Nq ~ A^^ the series in n in Eq. (!30ll 
with the first few terms, rapid convergence of the Taylor expansion of the exponential 
around the Gaussian form is assured by the following property of the matrix B{Q) 
defined in Eq. ( l32l) . Each term in the Taylor expansion of the right-hand side of Eq. ( l30l) 
around the Gaussian form {n = 2) generates a factor of the form 

J dQ Y[{B'''{n)) where J] fci = even = 2k and where all A;, > 3 . (36) 

i i 

That expression can also be written as 

N 
i fJ.=l 

In magnitude, each eigenvalue 6^(fi) is bounded by unity. It is, therefore, safe to expect 
that on averaging over the A^^g-dimensional unit sphere, each eigenvalue is of order 
l/y/N so that the expression in Eq. fl37j) is of order N~^VI{Nq). Here Vl{Nq) is the 
surface of the unit sphere in Nq dimensions. That shows that only few terms in the 
expansion are expected to contribute significantly even for Nq ~ A^^. The statement 
holds a fortiori for A^g <^ A^^. 

5. Proof 

To calculate the infiuence of the constraints in Eq. (12^ on the spectrum for A^ ^ 1, we 
use the approach developed in Ref. [8] based on the supersymmetry method [TT|, [12] . 
We only sketch the essential steps, using the definitions and notation of Refs. [12] and 
[8]. The average level density and all correlation functions are obtained with the help 
of a generating functional Z which is written as 

Z = j d^(exp{^^tL^/2QLi/2^|^^ ^3g^) 
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where 

G = H - E + M . (39) 

The average over the ensemble is denoted by angular brackets with an index H while for 
the trace we continue to use angular brackets without index as in Eq. ([2]). The symbol 
\E' stands for a supervector the dimension of which depends on the particular correlation 
function under study. The same is true of the matrices H (the Hamiltonian) , E (the 
energy), and M. The matrix M is of order 0{N~^) and contains energy differences, 
source terms and possible couplings to channels. Differentiation with respect to the 
source terms generates the particular correlation function of interest. 

The invariance of CGUE under unitary transformations implies that the integrand 
in Eq. ( |38ll depends upon and \E'^ only via the invariant form 

A^p = N-'lU^ ^ ■ (40) 

Here a amd f3 are matrix indices in superspace while fi runs over the N basis states of 
Hilbert space Ti. We introduce a supermatrix a with the same dimension and symmetry 
properties as A by writing Z as an integral over a delta function, 

Z = Jd^ Jda5{a-A){ exp {^^^V^^GV^^^^}) . (41) 

The delta function is replaced by its Fourier transform, and the multiple Gaussian 
integral over the supervector \E' is performed to yield 

Z = J da J dr exp{^A^ trg(r(T)}( exp { - ^tr trgln[G - Ar]})^ . (42) 

The remaining superintegrations over r and a are eventually done for N oo with 
the help of the saddle-point approximation. Prior to that step, the average over the 
ensemble is performed using Eq. (126|) . We first integrate over the unitary group. To 
leading order in N^^ we obtain 

(exp { — -tr trgln[G — Ar]})^^ = (exp { — -tr trglnD} 

X exp { - hi trg In [1 + j^tiD-^M]})^ . (43) 

Here D = x — E — Ar is diagonal in Hilbert space, and the angular brackets now stand 

for the remaining integration over the eigenvalues {x^}. Under inclusion of the terms 

which arise from W-p in Eq. (!26l) . the exponent of the integrand is given by 

111 
- -tr trglnD - -tr trgln [1 + — trD^^M] 

+ 2^1n|x^-x,|--^^xJ + lnFp(^) . (44) 

Following Ref. [8J, we perform the eigenvalue integration using the saddle-point 
approximation for ^ 1. In expression (jS]), all terms in the first line are at most of 
order while the first two terms in the second line are of order iV^. Omitting the terms 
in the first line but keeping In F-p(if) (depending on the value of Nq^ that function may 
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or may not be of order N"^), we put the derivatives of the resulting expression with 
respect to the equal to zero and obtain the N saddle-point equations 

2A2 ^1 2A2 1 dF-p{H) 



x„ = — y + , /i= l,...,iV . (45) 



Without the term InF-p(if) in expression (14^ . the saddle-point equations would have 
taken the standard GUE form 
2A^ ^ 1 

= ^ E ■ (46) 



To prepare for the treatment of Eq. fHS]) we recall how Eqs. fHB]) are solved in Ref. 
The variables x^/\ are replaced by a single dimensionless continuous variable, x^/\ e. 
The normalized level density p{£) with / de p{e) = 1 of the ensemble (which at this point 
is an unknown function) is introduced, and Eq. P6l) is written in the form 

e = 2P j de' . (47) 

Here P J stands for the principal-value integral. Eq. P71) has an electrostatic analogue 
and can be solved using the theory of analytic functions. The result is Wigner's semicirle 
law for p{e). The generating functional is subsequently taken at the saddle-point values 
for the Xfj_. All summations over in Z are, thus, replaced by integrations over e with 
p{e) as weight function. 

In applying that same method to Eqs. fH5|) we introduce the (yet unknown) 
normalized average level density p-p(e) of the constrained ensemble in the sum on the 
right-hand side of Eq. (115!) . We also have to implement the change of variables x^/ \ e 
in F'p{H) and its derivative. As for F-plH), this is done by replacing everywhere in 
Eq. ([30]) the term (1/A^)((i7/A)") by (e") = Jdee^'prie). The form of IV (i/) in Eq. (EEl) 
implies that pp{e) = p-p{—e) so that only terms with n even survive, and we obtain 



\27iNj 

X 



/ dt exp{Y^{l/{2nM~ir/N'--'){e'n{B'-m . (48) 

n>l 

This is a function of the unknown level density p-p with a rapidly converging Taylor 
expansion around the Gaussian term (n = 1). In the derivative of F-p, we substitute 
Xf_i/^ — ^ £ after differentiating with respect to x^. We obtain 

N^,-lrW) \2^) h^^^J ^ ^'^ 

X exp{j:{l/{2nm-ir/N'-'){e'-){B'-m . (49) 



The right-hand side of Eq. fl49l) is a polynomial of odd order in e with rapidly decreasing 
coefficients. 

As a result, the saddle-point equations (H5l) take the form of Eq. (147|) . with e 
replaced by an odd-order polynomial in e with rapidly decreasing coefficients. These 



Spectral fluctuation properties 



11 



coefficients depend on ppi^e); the solution of Eq. (H^ must, therefore, proceed iteratively, 
with the GUE average level density as a starting point for calculating (e"). In Ref. [1] 
it was shown perturbatively that p{E) and pv{E) differ by a term of order Nq/N'^. 
Therefore, we expect the two level densities to differ significantly when Nq ~ A^^. 

Returning to Eq. (jHj) we perform the integration over the variables a;^ by taking 
their values at the saddle points. That means, for instance, that we write 



This is the essential step: The summations over the eigenvalues disappear in all 
expressions in the integrand of Z. Each such summation is replaced by an energy 
integral involving the level density pvi^) of the constrained ensemble. This is the only 
place where the constraints show up in the calculation. From here on the calculation 
of the correlation functions for the CGUE and that for the GUE run completely in 
parallel [8] . It follows that all correlation functions of the CGUE have the same form as 
their GUE counterparts except that we have to replace the local average level spacing 
of the GUE by that of the CGUE. This proves our theorem. 

6. Discussion 

Our theorem holds in the limit ^ 1 and for Nq < Nq^'^. The average level density 
of the CGUE may differ from that of the GUE but all correlation functions have the 
same form for both ensembles. Although our result is perhaps expected, to the best 
of our knowledge this is the first time that GUE-type statistics has been analytically 
proved for a class of ensembles different from the GUE. Deviations of order 1 /N from the 
asymptotic form of the GUE statistics exist, of course, even for the pure GUE and are 
expected a fortiori for the CGUE. Our proof specifically applies to the case of unitary 
invariance. We believe, however, that a corresponding result holds also for the other 
symmetries. 

The proof of the theorem rests on the fact that in the limit A^ ^ 1 and for 
Nq < Nq^^, the constraining function Fp{H) is free of singularities. The proof holds 
independently of any specific properties of the constraining matrices Bg. What happens 
for A^ 1 but A'q > A'q'*? That seems to depend on specific properties of the 
constraints which determine the eigenvalues b^{s) and, thus, the convergence properties 
of the integrals over s. Therefore, generic statements about the spectral fluctuation 
properties of the CGUE probably cannot be made for Nq > A"q"*. 

One may speculate that with A^q increasing beyond the value Nq^^, the spectral 
fluctuations of the CGUE remain GUE-like until Np = N"^ — Nq is reduced to the value 
Np = N (where Hp = J2p hpBp may be a linear combination of A^ commuting matrices 
and, thus, integrable). But that speculation is surely incorrect. Indeed, random band 
matrices with a band width less than or of order -y/iV are known [161 [13 to possess 
localized eigenfunctions and a Poisson spectrum. For such matrices, the number Nq of 




(50) 
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constraints is at least of order A^^ — and, for ^ 1, obviously much larger than 

Nq but still much smaller than A^^ — A^. 

It is of interest to discuss the embedded random fc-body ensembles EGUE(A;) (see 
Ref. [3j and the review [1]) in the light of these considerations. The EGUE(k) models a 
Fermionic many-body system with fc-body interactions: m identical spinless Fermions 
occupy I degenerate single-particle states. The Hilbert space is spanned by A^ = 
Slater determinants. To construct the fc-body interaction operators, we denote by ajj 
and the creation and destruction operators for a Fermion in the single-particle state 
labelled /i with fi = 1, . . . , /. Let /ii, . . . , /i^ with 1 < /ii </ for all i denote a set 
of m non-equal integers, and analogously for z/i,...,z/m with 1 < z/j < / for all j. 
Then a general interaction operator has the form A{{fii}, {i^j}) = IliLi (^1,,^ YVJLi ^Uj- 
the Hilbert space of Slater determinants, the A^^ Hermitean operators A{{fii}, {i^j}) + 
^H{f^i}y {^j}) i[A{{fii}, {I'j}) — y4^({/ij}, {z^j})] play the very same role as do the 
matrices introduced in Section [2] in the Hilbert space 7i. From the general form 
of A, a /c-body operator is obtained by imposing the condition that a subset of m — k 
elements of the set {fii} is identically equal to a subset oim — k elements of the set {vj}- 
There are ('T") such fc-body operators. The EGUE(A;) is obtained by writing 

the Hamiltonian as a linear combination of all fc'-body operators with k' < k and with 
coefficients that are real random Gaussian-distributed variables. 

Among the EGUE(fc), the EGUE(2) has received particular attention because it 
mimics a Hamiltionian with two-body interactions, a form typical for Fermionic many- 
body systems like atoms or nuclei. One of the central questions (undecided so far) has 
been whether for A^ ^ 1 the spectral fluctuation properties of EGUE(2) are GUE- 
like. Numerical simulations [4] suggest that the answer is affirmative. However, these 
simulations are typically done for small values of / and m, with / around 12 and m 
around 4 or so. But for these values the number of one- plus two-body interaction 
terms is (\^) [4 x 8 + 6 x 28] = 200 (\^). The number of constraints is accordingly 

given by Nq = (^^^^ [495 — 200] = (^^^^ x 295. That figure is not much larger than 

A'q'* = (^^^ X 247, so that it is difficult to draw firm conclusions. It would be more 
informative to investigate numerically large values of / and m but that is prohibitively 
difficult. For Z ^ m ^ 1 we have N ^ l"^, and the number of independent two- 
body operators is approximately l^m^t^ . In other words, there are only m^/^ non-zero 
interaction matrix elements in every row and column of the matrix representation of 
the Hamiltonian for the EGUE(2), much fewer than for a banded random matrix where 
that number would be approximately = Z™/^. Put differently, the number Nq of 
constraints for the EGUE(2) is much bigger than it is for a banded random matrix. That 
fact suggests that mixing in the EGUE(2) is weaker than it is for a banded random 
matrix, and that EGUE(2) has Poissonian level statistics. On the other hand, in a 
banded random matrix it takes approximately different interaction matrix elements 
to connect two arbitrary states in Hilbert space. In the EGUE(2) that number is only 
m/2, i.e., less even than (1/2) In A^. This fact suggests that mixing of the states in Hilbert 
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space is much more efficient for the EGUE(2) than it is for a banded random matrix, 
and the question remains undecided. But the discussion suggests that for Nq > Nq^'^, 
the form of the constraints (and not just their sheer number) becomes important in 
determining the spectral fluctuation properties of CGUE. 

Another frequently used ensemble that simulates the nuclear many-body system is 
the two-body random ensemble (TBRE), see Refs. [5], [6] and the review [7]. Actually 
that ensemble is taken to be invariant under time reversal and, thus, has orthogonal 
rather than unitary symmetry. For simplicity we disregard this fact. The single-particle 
states belonging to a major shell of the nuclear shell model are occupied by a number of 
nucleons. The resulting Slater determinants are coupled to states with fixed total spin 
J and isospin T and are written as \JTfi). The running index n has a typical range R 
from several ten (J large) to several thousand or more (J small). Level statistics can 
be meaningfully discussed only for large R. It is assumed that the interaction between 
nucleons is of two-body type. Within a major shell, the number of independent two- 
body matrix elements Va is small (of order 10 or 10^) compared to the large values of 
R that are of interest. These matrix elements are taken to be uncorrelated Gaussian- 
distributed random variables. This defines the TBRE. For a given set of states | JT/x), 
the matrix representation of the Hamiltonian i^TB of the TBRE takes the form 



The matrices C^J(a) are fixed by the major shell and by the quantum numbers J and 
T under consideration but have some properties in common with matrices drawn from 
a canonical random-matrix ensemble. Again, it is of interest whether in the limit of 
R ^ 1 the TBRE generically obeys GUE (or GOE) level statistics. Numerical results 
and semi-analytical arguments both support such a hypothesis. Unfortunately, the 
matrices C^J(a) are not accessible analytically so far. Therefore, it does not even seem 
possible to determine the number Nq of constraints that would characterize the TBRE 
matrix ( 15T1) . and we cannot apply our results to that ensemble. 

The authors of Ref. [1| considered not only the CGUE but in addition also what 
they called "Deformed Gaussian Ensembles". Here the delta functions in Eq. (JT5l) are 
replaced by Gaussians, and the constraining function F-p{H) is everywhere regular. 
Following the arguments in Section [5] we conclude that the spectral fluctuations of the 
deformed ensembles coincide with those of the GUE. In other words, constraints affect 
the spectral fluctuation properties only if they constrain the relevant matrix elements to 
the value zero (and not to very small non-zero values). As a consequence, in the GUE 
the transition from GUE to Poisson level statistics is not a continuous process (where 
the level statistics would be smoothly deformed) but actually happens discontinuously. 
These statements apply on the "macroscopic" level where the values of the coefficients 
hq of the constraining matrices Bq are compared with those of the hp multiplying Bp. 
If, on the other hand, the hg are measured in units of the mean level spacing (i.e., 
in effect, on a scale l/\fN compared to the scale of the hp) then the transition from 
GUE to Poisson level statistics is expected to be smooth and to allow for intermediate 



{.H^:■B)^lv = 




(«)• 



(51) 
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forms of the level statistics. That expectation is supported by transitions between 
symmetries like the GOE GUE transition, and by many examples of partially chaotic 
systems that show intermediate level statistics. We have not attempted to introduce a 
correspondingly scaled parametrization for the deformed ensembles. Such a step would 
be meaningful only in the immediate vicinity of the transition point from GUE to Poisson 
level statistics. That point is not known analytically, however. 
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Appendix: Constraining matrices with non— zero traces 

It is convenient to relabel the indices q so that they run from 1 to Nq. We use 
Eq. (|T8l) and rotate the basis Bq Bq = Y.q' Oqq>Bqi with the help of an orthogonal 
transformation Oqg' such that Bi points in the direction of the unit matrix Ijy. Then, 

{Bq\Bq,) = 5qqi for all q' , 

{Bq) = for all g > 1 . (52) 

We apply the same orthogonal transformation to the variables tq so that tq ^ tq = 
J2q' Oqq'tqi and obtalu 

^^^^^ = (2^) y ndt.y ci[f/] exp(z^t,(5,|f//ff/t)) . (53) 

We write Bi as the sum of a traceless part and of a multiple of the unit matrix, 

B, = {B,-m'f) + m'f. (54) 

By construction, the traceless part is orthogonal to all the -Bg's with q> 1 and has norm 

{B, - m^\B, - (fii)^) = 1 - ^{B,f = a' . (55) 

Because of the first of Eqs. (I52l) . the eigenvalues bi^ with yU = 1, . . . , of i?i obey 
J2f_ibf^ = 1. We maximize (Bi) = X^^^ia* under that constraint and find that 

-Vn < (Bi) < Vn . (56) 

Therefore, 

< < 1 . (57) 
We define a as the positive root of and write Eq. (IMI) in the form 

5i = a5i + (5i)^ . (58) 
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Then, Bi has trace zero and norm one. We define = Bg for all g > 1 and have 

(Bg) = for all q and {Bg\Bgf) = 6qq> for all q, q' . (59) 

For a = the matrix Bi is a multiple of the unit matrix, and the integral over ii in 
Eq. (l53l) yields a multiple of the delta function for (H) while the remaining Nq — 1 
integrations over the ig with g > 1 are treated as in Section |H For a = 1 the matrix Bi 
is actually traceless; that case was treated in Section HI Therefore, we consider a only 
in the open interval 

< a < 1 . (60) 

We rewrite Eq. (!53l) by using the decomposition (!58|) . by rescaling aii ii, by 
introducing spherical polar coordinates {t, Q} in Nq dimensions, and by defining the 
matrix B{Q) as 

tB{n)=Y.igBg . (61) 

The function F-p^H) takes the form 

/ X'^ \Nq/2i f 

MH) = [^) -J dt t''^~'dQexpi^i^{t,n){B,){H)/{aN)) 

X Jd[U] exp{it{B{n)\UHU^)) . (62) 

Here ti{t,Q) stands for the rescaled old integration variable ti as expressed in terms of 
the new integration variables {t,Q}. The function F-p{H) diverges for (H) 0. The 
divergence is removed by multiplying F-p^H) with \{H / X)'^]^'^/'^ . As in the case of the 
substitution Fp — > F-p in Section [3], we expect that this step does not affect the spectral 
fiuctuations of the ensemble. The matrix B{Q) is traceless by construction, and the 
integration over the unitary group can be carried out as in Section |H From here on we 
proceed as in Section [51 
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